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Summary: A mathematical model is presented for the transport and precipitation of 
mineral in refilling osteons. One goal of this model was to explain calcification “halos,” in 
which the bone near the haversian canal is more highly mineralized than the more 
peripheral lamellae, which have been mineralizing longer. It was assumed that the precip- 
itation rate of mineral is proportional to the difference between the local concentration 
of calcium ions and an equilibrium concentration and that the transport of ions is by either 
diffusion or some other concentration gradient-dependent process. Transport of ions was 
assumed to be slowed by the accumulation of mineral in the matrix along the transport 
path. The model also mimics bone apposition, slowing of apposition during refilling, and 
mineralization lag time. It was found that simple diffusion cannot account for the transport 
of calcium ions into mineralizing bone, because the diffusion coefficient is two orders of 
magnitude too low. If a more rapid concentration gradient-driven means of transport exists, 
the model demonstrates that osteonal geometry and variable rate of refilling work together 
to produce calcification halos, as well as the primary and secondary calcification effect 
reported in the literature. 


In the last 50 years, the mineralization of bone has 
received much attention. Historically, many investi- 
gators have focused on the difficulty in the initiation 
of calcification caused by the metastability of inter- 
stitial fluid with respect to hydroxyapatite. Others 
have directed their attention to the structure of the 
completed mineral-organic complex. A foundation 
of basic concepts concerning the initial calcification 
processes and mineral-collagen structures has been 
built (2). 

Events occurring between the initiation of cal- 
cification and the completion of the fully calcified 
structure have received less attention. So-called pri- 
mary calcification takes only a few days, but the final 
30% or so of the mineral takes months to put in 
place (1,9,16). Little is known about the means by 
which calcium is transported from the bone surface 
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through layers of partially mineralized matrix to 
complete calcification of the still deeper lamellae. 
Diffusion is at least implicitly assumed to play a role 
(13), but it also is widely thought that cells actively 
transport ions during mineralization (16,19). 

Also, earlier concepts of calcium homeostasis re- 
cently have begun to give way to theories based on 
the exchange of bone mineral exclusive of remodel- 
ing (15). These theories cannot be fully elaborated 
until the transport, precipitation, and equilibrium of 
bone mineral are understood more completely. 

Observation of osteonal bone may provide in- 
sights regarding the transport of mineral during com- 
pletion of calcification. Microradiographs illustrate 
not only variations in mineral distribution between 
old and new osteons but also patterns of distribu- 
tion within the osteonal wall. Figure 1 shows that 
the innermost osteonal lamellae of newly formed 
osteons commonly are more radiopaque than lamel- 
lae farther out, which were laid down earlier and 
have been mineralizing longer. Something about the 
transport and precipitation of mineral tends to pro- 
duce a calcification “halo” adjacent to the canal 
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FIG. 1. Microradiograph of osleonal bone from rhesus monkey, showing variability in mineralization among osteons and the common 
occurrence of a “halo" of more heavily mineralized bone near the canal surface. The approximate field width is 1,000 pm Inset: 
Osteons just beginning to refill lack such a halo. 


(Fig. 2). Mineral density is high near the haversian 
canal, low in the middle lamellae, and increases again 
near the cement line. Weaver (22) established that 
the hardness of bone correlates with the halo; it is 
elevated near the haversian canal. 

I present here a mathematical model for the trans- 
port and precipitation of mineral in osteons; the 
model is intended to be a tool for achieving a better 
understanding of the middle ground between the ini- 
tiation of calcification and the completion of the cal- 
cified matrix. It addresses the following questions: 
(a) Can classic diffusion and precipitation equations 
serve as the basis for a finite-element model that 
converges on stable solutions and completes calcifi- 
cation in the observed period of time? (b) Does such 
a model produce calcification halos? (c) Does such 
a model exhibit primary and secondary mineraliza- 
tion behavior? (d) How is osteonal calcification af- 
fected by the geometry and dynamics of osteonal 
refilling, as opposed to crystal formation per sel 

METHODS 

Assumptions 

The present model does not deal directly with the 
complex and still unresolved issues surrounding the 
nucleation and growth of bone mineral crystals; it is 
concerned with the accumulation of mineral in the 


bone matrix once nucleation has occurred. Mineral- 
ization may involve the precipitation of a sequence 
of calcium phosphates having different solubilities 
(5,12). However, for the sake of the construction of 
an initial model, it was assumed that a single mineral 
(hydroxyapatite) is formed and that its rate ol pre- 



F1G. 2. Graph showing the result of scanning across the micro 
radiographic image of a mineralizing osteon. Radiographic den 
sity is plotted against scan position. The central depression is 
the haversian canal. There is a relatively high mineral content 
near the wall of the haversian canal and near the cement line 
The interior of the wall is less mineralized. Modified from 
Amprino and Engsirom ( t ) 
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cipitation is proportional to the difference between 
the local concentration of a single ion (calcium) and 
its equilibrium concentration (15). 

The model assumes that transport of ions is by 
either diffusion or another concentration gradient- 
driven mechanism. For example, transport of cal- 
cium ions via the osteocyte-canalicular network 
could be driven by the difference in the concentra- 
tion of this ion across the network. Thus, in what 
follows, the word diffusion should not be interpre- 
ted literally; indeed, the analysis will demonstrate 
that simple diffusion cannot be the primary trans- 
port mechanism. 

With this said, the model is based on Fick's Laws 
for diffusion, and it was assumed that transport 
of ions is controlled by the concentration gradient 
across, and the “diffusion” coefficient within, the ma- 
trix through which the ions must pass. The govern- 
ing equation is written in two-dimensional form so 
that curved osteonal bone-forming surfaces can be 
treated: 
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Here, C is the concentration of a specific ion species; 
D is an isotropic diffusion (or transport) coefficient 
for this ion; x and y are directions normal and paral- 
lel to the mineralizing surface, respectively; t is time; 
and P is the rate of precipitation of mineral ions from 
solution. For the sake of simplicity, it is assumed that 
Eq. 1 describes the transport and precipitation of 
calcium ions alone. 

The concentration is given as milligrams of free 
calcium per milliliter of interstitial fluid, and the rate 
of precipitation is given as milligrams of calcium 
precipitated per milliliter of interstitial fluid per sec- 
ond. To simplify matters, the rate of precipitation 
includes the nucleation of new crystals and growth 
of existing crystals. The rate is assumed to be pro- 
portional to the difference between the local con- 
centration of calcium ions (C) and an equilibrium 
concentration (C cq ): 

P = k cal (C - C eq ) n (2) 

where k ca , is a rate coefficient and n is a small number 
(12,15). This relationship is not limited to precipita- 
tion; it is assumed that if C < C eq , mineral crystals will 
dissolve. 

Calcification (I ) is defined as milligrams of cal- 
cium in hydroxyapatite per milliliter of bone. Typi- 
cally, 1 ml of cortical bone contains about 2 g of 
calcified matrix, of which 1.26 g is mineral, primarily 


hydroxyapatite (23). Since, by weight, 40% of hy- 
droxyapatite is calcium, the average calcium content 
of cortical bone, F*, is approximately 500 mg/ml. 

Microporosity is defined as the water-filled space 
amidst the collagen molecules and mineral crystals 
in the calcified matrix. The microporosity of osteoid 
contains the space to be filled by mineral, as well as 
the water that is bound to collagen or mineral, or 
both, and cannot be displaced. The volume fraction 
of the micropores is denoted by V w . 

The diffusion and precipitation of ions depend on 
the concentration of ions in the interstitial fluid 
within the micropores of the mineralizing matrix. 
Therefore, ionic concentration needs to be expressed 
per unit volume of interstitial fluid, not per unit vol- 
ume of bone. On the other hand, the rate of accumu- 
lation of crystalline mineral per cubic millimeter of 
bone clearly depends on the volume fraction of the 
micropores; if this volume fraction is small, there will 
be less interstitial fluid from which mineral ions can 
be precipitated into crystals, and the rate of calcifi- 
cation will be small relative to the rate of precipita- 
tion. Thus, the rate of accumulation of crystalline 
calcium, dT/dt, is assumed to be 

dF/dt = PV U (3) 

It is assumed that the volume fraction of mineral 
is zero in osteoid and that the organic fraction is the 
same in osteoid as in fully mineralized bone. As os- 
teoid mineralizes, mineral displaces interstitial fluid, 
reducing its volume fraction from about 66% to 
about 28% (23). This greatly diminishes the solute 
space through which mineral ions diffuse to reach 
growing hydroxyapatite crystals. The diffusion coef- 
ficient is assumed to diminish as mineral accumu- 
lates in the matrix, slowing transport. Neuman (13) 
hypothesized this phenomenon as follows. "As the 
crystals grow, they must grow closer and closer to- 
gether, the hydration shells must ultimately interact, 
and the diffusion of ions through the highly concen- 
trated jackets. . . becomes slower and slower. As the 
crystal faces approach each another, the restriction 
of charge and volume gradually will slow diffusion 
to a point where it is stopped almost completely." 

Thus, as calcification increases, it is assumed that 
the diffusion coefficient (D) decreases. A diffusion- 
calcification coefficient, k dc , is postulated such that 

D = D 0 exp(-k dc D (4) 

where D 0 is the diffusion coefficient for osteoid. 

Finally, the model incorporates two observations 
about osteonal refilling. First, as an osteon refills, the 
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area of the bone-forming surface (through which 
ions must pass to reach deeper lamellae) growls 
steadily smaller, "choking” the supply of mineral. 
Also, the innermost lamellae are laid down more 
slowly than the outer lamellae (8,10), giving them 
more time in close proximity to the source of ions. 

Model 

The model is based on the simultaneous applica- 
tion of the equations to an idealized bone-forming 
site (Fig. 3). Bone is being laid down by osteoblasts. 



CL MF OS 


FIG. 3. Schematic diagram showing a developing wall of new 
bone. Osteoblasts (dark ovals) are advancing the osteoid sur- 
face (OS) in the x direction with velocity M. Behind the osteoid 
(white region), a mineralization front (MF) follows at the same 
velocity. Bone formation began on a cement line (CL). The 
evenly spaced vertical lines show the relative x dimension of the 
elements used in the finite-element mesh. The completed osteo- 
nal wall contained 60 such elements. 

starting on a cement line. The bone behind the ce- 
ment line is assumed to be fully mineralized, but its 
mineral cannot move into the region of new bone. 
This is justified by the observation that, on microra- 
diographs (Fig. 1), a sharp line of demarcation is seen 
at the cement line when a partially calcified new 
osteon lies within a region of fully calcified bone. 
There is no evidence that mineral "leaches out” of 
the old matrix and crosses the cement line (1). 

The apposition rate of osteoid (M, in pm/day) is in 
the x direction. The surface of the osteoid in contact 
with the osteoblasts is termed the osteoid surface. At 
the mineralization lag time (T lil? , in days), after a 
layer of osteoid has been laid down, it begins to 
mineralize, forming a mineralization surface. The 
mineralization surface subsequently follows the os- 
teoblasts away from the cement line, moving to the 
right with velocity M. When bone formation ceases, 
the distance between the cement line and the osteoid 
surface is the wall thickness (X w .,„) of the completed 
formation site. 

Taken together, these equations constitute a non- 
linear system for which a closed-form solution would 


be very difficult. Therefore, a computer program in 
compiled BASIC (BVT0-5) by Bear and Verruijt (3) 
was modified to solve Eq. 1 while satisfying the other 
equations as well. Their program uses a central fi- 
nite-difference method to determine the transient 
behavior of a finite-element model for Eq. 1 They 
have described the details of the finite-element solu- 
tion and verification of its convergence to a closed- 
form solution (3). Their program was modified to 
simulate the continuous apposition of osteoid and 
other assumptions as follows. 

A mesh containing N elements, each X ctcm = 
X waM /N mm thick (x dimension), is defined (Fig. 4). 
To simulate refilling of a segment of an osteon, a 
wedge-shaped mesh is created by letting the y di- 
mension of each element linearly diminish so that 
the heights of the first and last elements are in pro- 
portion to the ratio of the radius of the cement line 
to the radius of the haversian canal. For other testing, 
a rectangular mesh is used. 

The mesh defines a completed wall of new bone. 
Initially, however, the model consists of only the left- 
most element, with its right edge nodes fixed at C 0 
to represent the presence of a vascular supply of 



FIG. 4. Drawn-to-scale illustrations of the meshes used in the 
finite-element model. The top, wedge-shaped mesh containing 
60 elements, was used for the osteonal model. Each element is 
1.333 pm wide. The heights of the first (farthest left) and last 
elements are 50 and 10 pm, respectively. The lower mesh, repre- 
senting formation on a flat surface, was the 60 element model 
used in the convergence test. Its overall dimensions were 80 x 
10 pm. Other meshes used in the test were of the same overall 
width and height but contained 30, 40. 50, and 70 elements of 
uniform width. 

calcium ions at the osteoid surface. Its left edge 
nodes are initially set equal to C n . but allowed to 
vary, and its diffusion coefficient is set equal to D„. 
This starting model is solved incrementally over the 
time required to lay down a single element (X t ., L . m /M). 
At the end of each time increment (§t), the rate ol 
precipitation (P) within the element is calculated 
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TABLE L Baseline values of model variables 


D () - 0. 10 mnr/sec 
K a . = 0.(X)1 ml/mg 

C u = 0.05 mg/ml 

C lN = 0.01 mg/ml 

K,. a i = 2.0 mmVmg-sec 
n = 2 

T* = 500 mg/ml 


d min = 3,300 mg/ml 
V or(t = 0.34 

Ti a( , = 10 days 
K fl |, = 0.02 day ‘ 


Diffusion coefficient for osteoid 

Diffusion-calcification 

coefficient 

Calcium concentration at free 
surface 

Equilibrium calcium 
concentration 

Crystal growth rale coefficient 

Crystal growth rale exponent 

Calcification of fully 
mineralized bone 

Density of bone mineral 

Volume fraction of organic 
phase 

Mineralization lag lime 

Nonlinear refilling coefficient 


from the mean of the element’s nodal concentrations 
by Eq. 2. The calcification of the element is calcu- 
lated as 

r = r + PV w 5t (5) 

where F is the calcification at the end of the pre- 
vious time period, PV w 8t is the amount of calcium 
precipitated per milliliter of bone during the current 
time increment, and V w = 1 - V org - V min . Here, V tirg is 
the volume fraction of the organic part of the bone, 
taken to be constant at 0.34 (23). V min is the mineral 
volume fraction, taken to be I7(m ca d min ), where m ca = 
0.40 mg of calcium per milligram of hydroxyapatite 
and d min = 3,300 mg/ml, or the specific gravity of 
hydroxyapatite (21). Then, the diffusion coefficient 
of the element is calculated from the calcification (T) 
by Eq. 4, and the model is solved again for the next 
time increment. 

When it is time to add another element to the 
model (X dcm /M = about 1 day), the mesh is rede- 
fined to include a second element; the right nodes of 
the new element are fixed at C 0 , and the left nodes 
(which had been the fixed right nodes of the original 
element) are made variable. The original element 
continues with the same initial values for the con- 
centration at its nodes and for the diffusion co- 
efficient and calcification (T) that it had at the end 
of the previous time period. (This is a mixed initial- 
value/boundary-value problem.) The new element is 
assigned D = D 0 and r = 0. This revised model then 
is solved through another X ctcm /M time interval. By 
continuation of this process, all N elements are ad- 
ded to the model. After the last clement has been 
added, the model is run an additional 200 days to 
approximate the completion of mineralization. 


To account for the mineralization lag time, ele- 
ments are not allowed to begin precipitating cal- 
cium until they have been in the model for T (ag days, 
creating an osteoid seam, where k ca , = P = V = 0 and 
D = D () . Precipitation was halted when V - I'*, de- 
fined as complete mineralization. 

The rate of apposition decreases during osteonal 
refilling (8,1 1) according to 

M = dR/dt = k fill R (6) 

where k fil[ is a refilling rate coefficient. To model this 
effect, R is assumed to be equal to the radius of the 
cement line minus the thickness of the wall added 
thus far, and the rate of apposition is recalculated for 
each new element. 

To determine an appropriate element size, a se- 
ries of rectangular (Fig. 4), constant rate of apposi- 
tion models of increasing resolution was tested with 
use of the model’s parameters (Table 1). Figure 5 
shows the results of this convergence test. The ef- 
fect of element size rapidly diminishes as N ap- 
proaches 70 elements. To restrict the solution time 
to about 1 hour (on a 486 personal computer oper- 
ated at 33 MHz), N = 60 was chosen for the analy- 
ses reported here. 

The speed and accuracy of the model also are 
affected by 8t. For accuracy, when a new element is 



POSITION, microns 

FIG. 5. Calcification plotted against position in the wall of 
newly formed bone for five models containing increasing num- 
bers of elements. This test of the effect of element size and 
number on the behavior of the model used a model configured 
for bone formation on a flat surface with a constant apposition 
rate (k ca , = 1.0 mmVmg-sec) and a mineralization lag time of 0 
days. 
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added, the first 5t should be <(X clcm ) : /2D 0 « 10 5 sec 
(3). Fortunately, 5t need not remain this small be- 
cause the rate of change within the model diminishes 
with time (3). Therefore, 5t is increased gradually, 
with the sequence (in seconds) of 0.0000 1, 0.00003, 
0.0(X) 1 , 0.0003, 0.001 , 0.003, and so on. Testing showed 
that this did not affect the results significantly. 

Selection of Parametric Values 

Values to be used for the variables in the model 
must be estimated. For lack of data, these baseline 
values often were chosen rather arbitrarily, but in 
each case they were tested by varying them para- 
metrically in the model until each gave results con- 
sistent with observed temporal and spatial patterns 
of mineralization. 

The baseline values for the mineralization lag time 
(10 days), radius of the cement line (R () , 0.100 mm), 
radius of the haversian canal (R h ,0.020 mm), refilling 
time (T fi ||, 80 days), and completed wall thickness 
(R c - R h , 0.080 mm) are typical of values in the liter- 
ature (11). The refilling coefficient was calculated 
from Eq. 6 as k fl „ = -ln(R h /R 0 )/T fm = 0.020 day 1 , with 
use of the baseline values. The equivalent mean rate 
of apposition was 0.001 mm/day. For convenience, 
complete calcification (F*), was set equal to the aver- 
age calcification (500 mg/ml) calculated earlier for 
cortical bone. 

The concentration of calcium in blood serum is 
about 2.48 m M or 0.10 mg/ml (7). Of this calcium, 
about 48% is free (not bound to proteins). The con- 
centration of calcium ions in interstitial fluid is 
thought to be somewhat more than half of that in 
blood plasma (20), but a higher proportion is free. 
Therefore, a reasonable estimate of the concentra- 
tion of calcium at the osteoid surface is 0.05 mg/ml. 

No data were found for the diffusion of calcium or 
phosphate ions through collagenous tissue or bone. 
The measurement of the diffusion of these ions 
through bone is confounded by the fact that when a 
concentration gradient is placed across the tissue, 
these ions are apt to precipitate within, or dissolve 
out of, the bone, as well as diffuse through it (14). 
Therefore, it is necessary to estimate the diffusion 
coefficient indirectly. 

The diffusion coefficient for strong electrolyte 
ions (such as calcium) in water at room temperature 
is slightly larger than 1 x 10^ cnr/sec (21). Roskey 
(6) determined the diffusion coefficient for calcium 
ions in denatured collagen gels to be 6.0 ± 0.5 x 10 * 
cm 2 /sec. On the basis of these limited observations, a 
reasonable baseline value for the diffusion coeffi- 


cient of calcium ions through osteoid (D„) would be 
5x10 6 cnr/sec = 0.0005 mnr/sec. However, when 
this value was tried in the model, it was found to be 
much too low to produce observed levels of calcifi- 
cation. After testing of various values, D 0 = 0.10 
mm 2 /sec was selected as a baseline value. Also, lack- 
ing experimental data, 0.001 ml/mg was assigned as 
a baseline value for the diffusion-calcification coef- 
ficient (k dc ) after testing of various values in the 
model. 

To learn how well Eq. 3 applies to the precipita- 
tion of hydroxyapatite, an experiment by Boskey 
and Posner (4), and reported on by Wuthier (24), was 
analyzed. In this experiment, hydroxyapatite crys- 
tals were grown in a solution without replenishment 
of the calcium (or phosphate) ions. An equation de- 
scribing this situation was derived with use of Eq. 3 
and found to fit the experimental data very well 
when n = 2. The result also showed that the equi- 
librium concentration (C cq ) was 0.025 mg/ml if the 
precipitating “seed'' was a phospholipid, and 0.006 
mg/ml if it was hydroxyapatite crystals. Correspond- 
ing values of the rate coefficient (k cal ) were assumed 
inappropriate for the model because nucleation sites 
were probably much more sparse in the experiment. 
After different values of the rate coefficient and the 
equilibrium concentration were tried in the model, 
baseline values of 2.0 ml/mg-sec and 0.01 mg/ml, re- 
spectively, were selected. 



WALL POSITION, micro 

FIG. 6. Calcification profiles (F compared with position in the 
osteonal wall) are shown 30, 80, 123, 168. and 280 days alter the 
start of refilling at a variable rate with baseline parametric values. 
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TIME, days 

FIG. 7. Graph showing temporal changes in calcification at 
three representative locations in an ostconal wall: near the 
cement line (5 pm), in the middle (39 pm), and near the haver- 
sian canal (74 pm). The first two curves exhibit primary (steep 
slope) followed by secondary (reduced slope) mineralization 
behavior. 

RESULTS 

Extensive testing of the model showed that when 
its parameters were kept within the vicinity of the 
baseline values, it was stable and insensitive to re- 
finements of element size and time increments. How- 
ever, as previously mentioned, the model does not 
support the concept of transport by simple diffusion 
because, no matter how much time is allowed or how 
other parameters are adjusted, the values for diffu- 
sion of calcium ions in collagen gels did not produce 
physiological levels of calcification. As 0.03 mm 2 /sec 
was exceeded, the model began to approximate 
physiologic results, and values in the vicinity of the 
baseline of 0.10 mm 2 /sec produced nearly complete 
calcification at 280 days. 

Early in the refilling process (30 days), there was 
no indication of a calcification halo (Fig. 6). This 
concurs with the appearance of the osteon begin- 
ning to refill in Fig. 1. Before refilling ended, how- 
ever, a halo emerged, and in the newly completed 
osteon (80 days), the profile was quite similar to 
that in Fig. 2. There was a peak of calcification near 
the canal surface, a valley in the central region of 
the wall, and a return to high levels near the ce- 
ment line. As calcification continued after the cessa- 
tion of bone formation, the region nearest the canal 
surface continued to mineralize faster than deeper 
regions (145 days), and a state of ''complete miner- 


alization” slowly worked its way toward the cement 
line. 

The model also exhibited primary and secondary 
calcification (Fig. 7). In the lamellae near the cement 
line and in the middle of the wall, there was an initial 
rapid calcification followed by a much slower rate of 
mineral accumulation. This profile is similar to that 
of primary and secondary calcification as described 
by Parfitt (16). However, the model predicted that 
lamellae laid down near the end of osteonal refilling 
rapidly reach full calcification. 

The comparison of the behavior of “flat” and os- 
teonal models (data not shown), with and without 
variable rates of apposition, showed that calcifica- 
tion halos are caused primarily by the variable refill- 
ing rate and secondarily by the bottle-neck effect. 
Also, experimentation with the model showed that 
achievement of physiologic calcification levels de- 
pends on parameters being within a fairly narrow 
range. For example, varying the diffusion-calcifica- 
tion coefficient (k dc ) had little effect while the osteon 
was refilling, but values less or greater than the base- 
line value retarded completion of mineralization; in- 
creasing the rate coefficient (k ca) ) over the baseline 
value enhanced calcification during refilling but later 
resulted in reduced mid-wall calcification. 

DISCUSSION 

The goal of this work was to answer several ques- 
tions concerning a mathematical model of mineral- 
ization based on classic equations for diffusion and 
precipitation. Does such a model exhibit stable and 
convergent behavior and produce observed levels of 
calcification in physiologic periods of time? The an- 
swer to the first part of this question is yes. The 
answer to the second part also is yes, but only if the 
rate of ion transport is much greater than that pro- 
vided by simple diffusion. Does the model produce 
calcification halos and exhibit primary and second- 
ary mineralization rates? The answers to both these 
questions are definitely yes. Finally, is osteonal min- 
eralization affected by osteonal geometry and refill- 
ing dynamics? Indeed, these factors were found to 
be primarily responsible for halo formation. 

Schmidt, in 1933, as cited by Amprino and Eng- 
strom (1), attributed the greater degree of calcifica- 
tion of the inner osteonal lamellae to their slower 
rate of apposition, but the latter authors believed 
that this was immaterial; they postulated that the 
inner lamellae were more calcified because they 
were closer to the vascular supply, so that “gradually 
the content of calcium salts becomes greater in the 
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inner part of the osteon” (1). The model suggests that 
each view contains an element of the truth: inner 
lamellae become more highly mineralized because 
they are close to the blood supply for a longer time 
during refilling and thereafter. 

One of the most important results was the quanti- 
tative demonstration that the mineralization of bone 
cannot be sustained on the basis of simple diffusion 
as the primary means of transport of mineral ions. 
The model supports alternative hypotheses that cells 
or matrix structures augment movement of ions if 
such transport is concentration gradient-dependent. 
It is, for example, reasonable to postulate that osteo- 
blasts and osteocytes transport mineral ions through 
their canalicular processes according to the differ- 
ence in ionic concentration between two points in 
this network. 

The success of the model in predicting certain ba- 
sic observations constitutes an obvious strength. At 
a deeper level, this model, like all mathematical 
models, is valuable because it enables testing of re- 
lated hypotheses for internal consistency as a system. 
There obviously are many limitations to the model. 
It provides a numerical rather than a closed-form 
solution, and although care has been taken to ensure 
that the solutions presented here are not signifi- 
cantly affected by the stability and accuracy prob- 
lems that can affect such models, future versions also 
must be carefully checked for such errors. It is two- 
dimensional in space and one-dimensional in terms 
of the mineral ions that it considers. It also has not 
incorporated the hypothesized existence of precur- 
sors to hydroxyapatite during calcification. 

Although the lack of experimental data defining 
parametric values is another important limitation of 
the model, experimentation with the model showed 
that many of the parameters had to remain within 
rather narrow ranges to give physiologic results. This 
provides the basis for two kinds of important hy- 
potheses: hypotheses for experiments designed to 
determine values for parameters (e.g., does the rate 
coefficient [k ca ,] fall within the range suggested by 
the model?) and hypotheses explaining pathologic 
conditions (e.g., are certain osteomalacias marked by 
abnormal values for the rate coefficient [k ca ,]?). With 
respect to the latter, it has been shown that alumi- 
num reduces the growth rate of hydroxyapatite crys- 
tals (18). The effects of this on the calcification 
profiles of osteons could be modeled and the result- 
ing calcification profiles compared with microradio- 
graphs of osteons affected by aluminum-induced 
osteomalacia. 


It should be noted that some combinations of 
parametric values caused the model to become quite 
dynamic. Mineral precipitated at one point was dis- 
solved, moved to another region, and precipitated 
again. Neuman (13) stressed that there is a continu- 
ous exchange of ions back and forth from the solu- 
tion through the hydration shell to the crystal. The 
model supports this view that dissolution, exchange, 
and crystal growth are all part of a continuum, 
and relatively small changes in local conditions can 
move the system from one state to another. This, in 
turn, is consistent with recent theories of calcium 
homeostasis (17), which postulate the movement of 
mineral across inactive surfaces, with lining cells 
and osteocytes controlling equilibrium. 

It must be emphasized that this is a preliminary 
model that should be extended to produce more sig- 
nificant results relevant to current questions about 
the calcification of bone. Future improvements could 
include the transport of mineral ions through the 
osteoblast-canaliculi-osteocyte network, followed by 
diffusion into the matrix; development of a multiple- 
ion version based on calcium-phosphate solubility 
products; and the incorporation of precursor mineral 
phases. 
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